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We revisit the intersubband polaritonics — the branch of mesoscopic physics having a huge 
potential for optoelectronic applications in the infrared and terahertz domains. We show that 
contrary to the general opinion the Coulomb interactions play crucial role in the processes of light- 
matter coupling in the considered systems. We demonstrate that electron-electron and electron-hole 
interactions radically change the nature of the elementary excitations in these systems. We argue 
that intersubband polaritons represent the result of the coupling of a photonic mode with collective 
excitations, and not non-interacting electron-hole pairs as it was supposed in the previous works on 
the subject. 



I. INTRODUCTION 

Intersubband optical transitions in semiconductor 
quantum wells (QWs) are widely used in a variety of 
modern optoelectronic devices operating in a broad wave- 
length range spanning from mid- infrared to terahertz 
Their practical realization needs high radiative quantum 
efficiency. In this context, the implementation of the con- 
cepts of strong light-matter coupling is a promising tool 
to improve the functionality of the devices as compare to 
those operating in the weak-coupling regime.^^ 

The achievement of the strong coupling is possible if 
the absorbing media is placed inside a photonic cavity 
and coherent light-matter coupling overcomes the dissi- 
pative processes in the system. For intersubband tran- 
sitions the experimental realization of strong coupling 
regime was for the first time reported in the pioneering 
work of D. Dini et al^ The elementary excitations in this 
case have hybrid, half-light half-matter nature and are 
called intersubband polaritons. They have a number of 
peculiarities distinguishing them from conventional cav- 
ity polaritons formed by interband excitons. First, they 
are formed only in TM polarization, as optical selection 
rules prohibit the absorption of the TE mode in inter- 
subband transitions. Second, the strength of the coupling 
can be an important fraction of the photon energy, which 
makes possible the transition to so-called ultrastrong cou- 
pling regime^iS. The light-matter coupling constant and 
resulting Rabi splitting depend on the geometry of the 
QW and photonic cavity and the electronic density in 
the lowest energy subband^^, which opens a possibility to 
tune this parameter by application of the external gate 
voltage. 

The broad variety of the applications of intersubband 
polaritonics makes important the understanding of the 
nature of intersubband polaritons. The question is: what 
kind of the elementary excitation in a QW is coupled with 
a photonic mode and participate in the formation of the 
polariton doublet? The former can be devided into two 
categories: single-particle excitations (SPE) and collec- 
tive excitations, appearing from electron-electron inter- 
actions and absent in the non-interacting system. The 



earlier works devoted to theoretical description of the 
intersubband polaritons neglected Coulomb interactions 
completelyiSrJ^ and the main qualitative conclusion was 
that the formation of the polaritons is a result of the 
coupling between non-interacting electron-hole pairs and 
a cavity mode. Moreover, it was claimed that bosoniza- 
tion approach is valid for the description of unbounded 
fermion pairs.— 

The opinion that Coulomb effects play no substantial 
role in the intersubband polaritonics seems controversial, 
as their important role in photoabsorption of individual 
QWs (in the absence of a photonic cavity) is an estab- 
lished fact, studied extensively from both experimental 
and theoretical points of view.—"— Electron-electron in- 
teractions lead to the appearance of the collective exci- 
tation modes such as intersubband plasmon (ISP) which 
under certain conditions can give a dominant impact to 
the optical responseJ^"— 

The only work considering the role of many-electron 
interactions for intersubband transition in the microcav- 
ity known to us is a paper by M. Pereira3^. However, the 
influence of interactions on quasiparticlc spectrum was 
not the main focus of the paper. The qualitative aspects 
of strong light-matter interaction of collective intersub- 
band modes with microcavity photons remained unin- 
vestigated and the role of the many-body corrections of 
various types thus remained unclear. 

In the current paper we bridge this evident gap. 
We propose a semi-analytical way of the description of 
the coupling between intersubband excitations (single- 
particle or collective) to a cavity mode in terms of Feyn- 
man diagrams corresponding to the different physical 
processes in the system. This makes our calculations 
transparent and allows a simple qualitative interpreta- 
tion of the role of many body interactions in the consid- 
ered system. Moreover, since we sum up infinite series 
of the diagrams, our treatment is non-perturbative and 
includes all orders of the interaction. As well, it treats 
the resonant and anti-resonant terms in light-matter cou- 
pling Hamiltonian on equal footing and thus allows the 
description of the phenomenon of ultra-strong coupling 
as well. We show that Coulomb interactions play cru- 
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cial role in the intersubband polaritonics, especially for 
high electron concentrations necessary for the achieve- 
ment of strong coupling regime. We claim that in this 
case the intersubband polariton is formed due to the cou- 
pling of the collective excitation known as intersubband 
plasmon (ISP) with the cavity mode. In the opposite case 
of small electron concentrations we show that cxcitonic 
corrections become crucial. 



II. THE MODEL 

We consider a system with GaAs/AlGaAs quantum 
well (QW) embedded into microcavity in the configura- 
tion usually used for intersubband transitions with TM 
polarized light (Fig. [Ha,b)). As our goal here is to 
present qualitative results and not to perform a fit of any 
experimental data, we do not consider the case where the 
photonic cavity contains several QWs which is often used 
in order to obtain larger values of Rabi splitting and nec- 
essary for the achievement of the ultra-strong couplingj^ 
Our method, however, can be easily generalized for this 
configuration as well. We choose a doping of a single QW 
in such a way that only the lowest subband is filled by the 
electrons at T = while upper subbands remain empty 
(Fig. [He)). Later on we concentrate on T = case only. 

When quantum well is embedded into a microcavity, 
the photons interact continuously with electrons mov- 
ing them from fundamental subband to the first excited 
subband, thus creating electrons and holes which can in- 
teract with each other. This electron-hole pairs can then 
again disappear, re-emitting a photon. In diagrammatic 
language such process can be described by a polarization 
bubble n (Fig. The photon in a cavity thus becomes 
"dressed" by such bubbles, and its Green function G can 
be described as a sum of the terms containing one, two, 
three etc bubbles as it is shown in Fig. [2{a). The sum- 
mation up to the infinite order gives the Dyson equation, 
whose solution yields 



TM-polarized light 



G = 



Gq 



(1) 



where 11 denotes to the full polarization operator de- 
scribed by the intersubband bubble containing all pos- 
sible Coulomb interactions. In this expression Gq is a 
bare photon Green function. 



Go(w,(?) = 



2hwQ{q) 



ri^uj^ - h^Ljliq) + 2irujQ{q)' 
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where Wo('?) describes the cavity mode dispersion and F 
is the broadening of the photonic mode due to the finite 
lifetime (taken to be « lOps) and g is a matrix element 
of electron-photon interaction which reads^ 
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FIG. 1: (Color online) Geometry of the system, (a), GaAs 
quantum well (QW) placed into the microcavity created by 
distributed Bragg reflectors or air. Using the total internal 
reflection of light, TM-polarized beam travels through the 
structure exciting the excitations between upper and lower 
subbands. The current figure shows the creation of inter- 
subband plasmon - charge-density excitation, (b), Sketch of 
the QW of width L = 12.8 nm with plotted wave functions 
for fundamental and upper subbands. The separation energy 
between levels is (A = 100 meV). (c), Dispersions of elec- 
trons for two subbands with Fermi level energy for electron 
concentration riei = 1.25 x 10^^ cm~^ . The single-particle 
excitations are schematically described as blue arrows. The 
kp label denotes to the Fermi wave vector. 



where Lcav is cavity length, A is separation energy be- 
tween levels, Co and e are vacuum permittivity and rela- 
tive material dielectric constant, respectively, dio stands 
for the dipole matrix element of the transition and A is 
an area of the sample. 

The dispersion of elementary excitation in such system 
is determined by the poles of G{uj,q) and can be found 
by solving a transcendental equation 



1 



fG^Ii = 0. 



(4) 



The coefficient of the photoabsorption, a{q,uj), is pro- 
portional to the imaginary part of the full polarization 
operator of intersubband polariton accounting for multi- 
ple re-emissions and re-absorptions of the cavity photon 
npo;((7,w): 



(5) 



The diagrammatic representation of the equation for 
Iipoi{qTLo) is shown in Fig. [Hb) and yields 



n 



n 



pol 



1 - .g^GoH 



(6) 



3 



(a) 



cavity photon 



(b) 



1 

ripoi^ 






bare photon 

1 


1 





intersubband 
excitation 



teractions is straightforward and givei 



,18 







1 









(c) 



quasiparticle bare eiectron 





(d) x-hf; 



seif-energy Hartree Foci< 

FIG. 2: Diagrammatic representation of intersubband polari- 
tons. (a), The Dyson equation for microcavity photon inter- 
acting with intersubband quasiparticles 11 leading to the for- 
mation of polariton. The label g corresponds to quasiparticle- 
photon interaction constant, indices 1 and denote to excited 
and fundamental subband, respectively, (b), Series for ab- 
sorption by intersubband quasiparticle coupled to the cavity 
in the general form with accounting all many-body effects, 
(c), The Dyson equation for self-energy corrections of the 
Green function of the electron, (d), The Hartree- Fock self- 
energy which consists of direct (Hartree) and exchange (Fock) 
diagrams. 



As it is seen from the Eqs. (|4])-(l6]), all properties of 
intersubband polaritons can be determined if the expres- 
sion of the polarization operator of an individual QW 
n(g, uj) accounting for Coulomb interactions is known. In 
general, the calculation of this quantity is a tricky task 
which can be performed only in some particular cases, 
which we are now going to consider. 



III. NON-INTERACTING CASE AND 
HARTREE-FOCK APPROXIMATION 

This case has a methodological interest and represents 
a test for our approach, allowing to compare the results 
it gives with those obtained earlier in the Refs. [Tol - [l3j 
by using the bosonisation scheme. For non-interacting 
particles the calculation of the polarization operator Hq 
represented by a single bubble without any Coulomb in- 
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where E-^ 

are dispersions of electrons in the excited and fundamen- 
tal subbands, respectively, and rtok is the Fermi distribu- 
tion in the fundamental subband which can be replaced 
by a step-Hke function at T = 0. The quantity 7 rep- 
resents a non-radiative broadening of SPE. This integral 
can be calculated analytically (see Appendix A) . The use 
of Eqs. (I4])-(|6]) allows the determination of the disper- 
sions of the intersubband polaritons and photoabsorption 
of the system. It is instructive to consider the case 7 — > 
and assume that the transferred momentum of photon q 
is small as compare to the Fermi momentum of the elec- 
tron gas kp. In this case, the equation for the energies 
of the polariton modes Q reads 

{hjjj — huQ + ir){f}aj + hujo — iT){hLij — A) = 2ujQnAg'^{q). 

(8) 

If we are outside the ultrastrong coupling regime (wo 3> 
nAg'^{q)) this reduces to 



(hjjj — hujo + ir){haj — A) = nAg^{q), 



(9) 



which is nothing but the equation for two coupled har- 
monic oscillators corresponding to a cavity mode and 
single-particle excitations in the QW. The dispersions of 
intersubband polaritons deduced from this equation co- 
incide with those obtained in the earlier works 

One should note that in the mean field (Hartree-Fock) 
approximation the electron-electron corrections can be 
easily introduced into consideration without substantial 
modification of the formalism. No new collective exci- 
tations appear in this approach and the results remain 
qualitatively the same as for the non-interacting case. 
We thus consider both situations in the same section. 

In diagrammatic representation the Hartree-Fock ap- 
proximation results into renormalization of electron 
Green functions which can be described by the Dyson 
equation shown in Fig. HJc). Only the first order dia- 
grams corresponding to the direct Hartree term and ex- 
change Fock term (Fig. [Hd)) are retained in the expres- 
sion for the self-energy in this approach. The Hartree 
term diverges in the limit A ^ 00 but is compen- 
sated by the interaction of the electrons with positive 
backgrounds^ and the Fock exchange self-energy correc- 
tion can be written in the form 



ki 



ViQOt{\k - ki|)nki 



(10) 



where i = 0,1 and the index corresponds to the funda- 
mental subband and index 1 to the first excited subband, 
denotes the Fermi distribution in the fundamental 
band and sign " — " shows that the exchange interaction 
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decreases the energy. Vijki{q) denotes a matrix element 
of the Coulomb interaction which reads 



2eeoAq 



dzdz'U^)cp^{z')dpk{z')Uz)e-''\' 



(11) 

with indices i,j,k and / corresponding to the initial and 
final subbands from which particles interact and 4'i{z) 
being the envelope wave functions in the direction of the 
structure growth axisJ^i^i 

The calculation of the polarization operator accounting 
for the Hartree-Fock corrections can be done by substi- 
tuting the energies of bare electrons in the fundamental 
and first subbands by their rcnormalized values, calcu- 
lated as 



i;W(k) = i;«(k) 



(k) 



(12) 



This rcnormalization has the following consequences. 
First, the correction for populated fundamental subband 
is greater then for the empty first subband leading to the 
rcnormalization of the transition energy A. Accounting 
for the negative sign of exchange correction, one con- 
cludes that the effective gap A is increased. Second, in 
general the self-energy is a function of momentum which 
leads to the non-parabolicity of the renormalized disper- 
sions and contributes to the broadening of the absorption 
line. For intersubband polaritons the first effect shifts the 
anticrossing point in the region of larger momenta (sev- 
eral meV in the geometry we consider) and second leads 
to the slight decrease of the observed Rabi splitting. 



IV. INTERSUBBAND PLASMON-POLARITON 
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In general, electron-electron interactions can not be 
neglected and 11 ^ Ho. Their account is a compli- 
cated task. However, for high enough electron con- 
centrations (n ~ lO^^CTO"^) the polarization operator 
can be estimated using Random Phase Approximation 
(RPA) whose diagrammatic representation is shown in 
Fig. |3l^c).— — In this regime the system demonstrates 
the appearance of intersubband plasmon (ISP) — charge- 
density collective excitation arising from intersubband 
transitions. 

The summation of the infinite series of the diagrams 
represented in Fig. [3Kc) allows us to obtain a compact 
expression for the polarization operator corresponding to 
ISP: 



n 



ISP = 



1 — ^loiollo 



(13) 



The dispersion of intersubband plasmon is plotted in 
Fig. (nija). The absorption by intersubband plasmon is 
given by imaginary part of polarization 11 (Fig. [3l^b)). 
One sees that the system still has an absorption peak 
corresponding to the single-particle excitations. How- 
ever, another peak corresponding to the absorption of 



^ ISP dominates over SPE in 
(b), Absorption spectrum of 
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FIG. 3: (Color online) Intersubband plasmon. (a), The dis- 
persion of intersubband plasmon (ISP, red line) and single- 
particle excitations spectrum (SPE, violet parabolas) plotted 
for the long range of momentum. For the considered concen- 
tration Uei — 1.25 X 10^^ cm 
the small wave vectors range 

intersubband quasiparticles plotted for momentum q = 10' 
m"^ . While the intersubband plasmon peak is sharp [5 — 0.05 
meV) and high, the SPE continuum is much broader [5 — 0.77 
meV) and has smaller oscillator strength, (c), Random Phase 
Approximation (RPA) series for intersubband electron-hole 
pairs interacting by direct Coulomb matrix element Vioio- 
These diagrams correspond to the formation of intersubband 
plasmon. 



ISP appears. This peak is blueshifted by a value of de- 
polarization shift nVioio{q). It is much more intensive 
and narrow that the peak corresponding to SPE. There- 
fore, it is natural to suppose that after placing of the QW 
in a photonic cavity this peak will give main contribution 
to the formation of intersubband polariton. 

This conclusion is supported by the results of calcula- 
tions shown in Fig. [H The dispersion of the elementary 
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FIG. 4: (Color online) Intersubband plasmon polariton. (a), 
Density plot of the intersubband plasmon polariton spectrum 
showing the dispersion of excitations and absorption (color 
intensity) in the system. The detuning of cavity mode is 10 
meV. (b), Absorption spectrum of the system plotted for wave 
vector q = 1.25 X 10® where the anticrossing point for 

ISP-photon exists. The two peaks corresponding to upper and 
lower plasmon polaritons are clearly observed, while single- 
particle excitations are suppressed, (c), Dispersions of inter- 
subband plasmon polariton modes (red and violet lines) in 
strong coupling regime with corresponding anti-crossing and 
Rabi frequency Vr = 2.4 meV. The SPE dispersion (pink 
line) is weakly coupled to the cavity mode. 



excitations of the hybrid QW-cavity system was calcu- 
lated using Eq. ([4]), where the Hartree-Fock corrections 
were accounted for in XIq. The corresponding spectrum 
of intersubband excitations in the semiconductor micro- 
cavity contains three branches (Fig. |l{c)). Two of them 
corresponding to the cavity photons and ISP reveal an- 
ticrossing and give birth to the intersubband polariton 
modes, which in this case can be called more correctly in- 
tersubband plasmon-polaritons. The third mode denotes 
to the single-particle excitations. The corresponding dis- 
persion line crosses the dispersion of the cavity photon, 
which means that for SPE the weak coupling regime is 
realized. In the absorption spectra shown in Fig. IH^b) 
three peaks appear. Two of them corresponding to inter- 
subband plasmon-polaritons are very pronounced and the 
third one corresponding to coupling with single-particle 
excitations is very weak but still observable (Fig. lU^a)). 
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FIG. 5: Diagrammatic representation of intersubband exci- 
ton. (a), The ladder series for the intersubband excitation 
leading to the creation of intersubband exciton. The wavy line 
corresponds to Viooi Coulomb interaction and double wavy 
line denotes to effective interaction between particles, (b). 
Two-particle integral equation for effective interaction, q is 
the transferred momentum, (c). The general Dyson equation 
which describes the situation of mixed excitonic and plas- 
monic effects leading to the formation of mixed collective 
modes. 



ing. In this case plasmonic corrections can be neglected 
and excitonic effects corresponding to the interaction be- 
tween the photoexcited electron in the first subband with 
a hole in fundamental subband become dominant. In di- 
agrammatic language they can be described by ladder di- 
agrams shown in Fig. ^a)^^^ Summation of the ladder 
diagrams up to an infinite order gives birth to the forma- 
tion of the attraction between the electron and hole, and 
can give a peak i photoabsorption lying below the con- 
tinuum of SPE (contrary to the the peak corresponding 
to ISP). 

Mathematically, the dispersion of the excitonic mode 
can be deduced from the integral Bcthe-Salpetcr equa- 
tion for the effective electron-hole interaction, whose dia- 
grammatic representation is given in Fig. and which 
reads2^ 



V. EXCITONIC EFFECTS 

In the previous section we investigated the case of high 
electron concentration in QW where plasmonic effects 
dominate and RPA can be successfully used. However, 
the opposite limit of low concentrations is also interest- 



iy(k,ki,L^,q) = -Viooi(k-ki)- (14) 
- J dk2Viooi{'k - k2)no(a;,k2,q)W^(k2,ki,a;,q). 

The solution of the Bethe-Salpeter equation is a compli- 
cated task which goes beyond the scope of the present 
paper. Here we restrict ourselves by a simplified treat- 
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FIG. 6: (Color online) Intersubband exciton polariton. (a), 
The dispersion of intersubband exciton (green line) which 
lies below the single-particle excitations continuum (violet 
parabolas). The concentration of electrons in QW is n^i = 
10^^ cm~^ . (b), Density plot of intersubband exciton coupled 
to the photonic mode. The intensity describes the absorp- 
tion in the system, (c), Dispersions of intersubband exciton 
polariton modes (green and violet lines) in strong coupling 
regime [Vr = 0.5 meV). The SPE dispersion (pink line) is 
weakly coupled to a cavity mode. 



mcnt of the screened 2D Coulomb potcntia' 



Viooi(q) 



24.27 



2eeoA(|q| 



(15) 



where k = „ ""^ is the screening wave vector. In 2D 
systems it does not depend on electron density and for 
the structure we consider is equal to k = 2.4 x 10*m~^. 
Then, since we are working with zero temperatures and 
small concentrations {kp < n), the interaction between 
particles is slowly varying function of q and can be ap- 
proximated by a constant value Vq w /2eeoAK. In this 
limit the equation for the polarization operator corre- 
sponding to the intersubband exciton can be found ana- 
lytically asi^ 



no 



1 + VoUq 



(16) 



where the negative sign of Coulomb interaction corre- 
sponding to the electron-hole attraction is accounted for 
in the denominator. The calculated dispersion of inter- 
subband exciton is shown in Fig. [D^a) for the QW with 
doping n^i ~ 10^^ cm^^, where assumption fc^^ < k is 
appropriate. 

The dispersion of the elementary excitations of the QW 
coupled to a cavity mode accounting for the excitonic ef- 
fects are shown in Fig. ^c). Similarly to the case of the 



high concentrations with the strong ISP-photon coupling, 
one sees three dispersion branches. Two of them corre- 
sponding to the excitonic and photonic modes reveal an- 
ticrossing and form the intersubband exciton polaritons. 
The third one corresponding to single-particle excitations 
remains in a weak coupling regime. In the photoabsorp- 
tion spectrum shown in Fig. IHJb) three peaks of different 
intensities are observed. The excitonic effects lead to the 
redistribution of the oscillator strength, which becomes 
small for SPE transitions and corresponding peak is con- 
sequently very weak. In this case the Rabi splitting is 
reduced due to the concentration dependence of interac- 
tion constant and is about 0.5 meV. 

One should also comment on the case of the interme- 
diate densities, where strictly speaking neither RPA nor 
ladder approximation can be applied. In some worka^S 
RPA and ladder corrections were accounted simultane- 
ously to describe the Coulomb correlations in this regime. 
In the diagrammatic representation this results into the 
equation for the polarization operator shown in Fig. [SJc) 
which leads to the appearance in the system of some 
mixed collective exciton-plasmon modes. As depolariza- 
tion and excitonic shifts have opposite signs, it may be 
possible that there exists a regime when they almost com- 
pensate each other and the energy of the collective mode 
is close to the energy of the single-particle transition. 
The consideration of this case, however, goes beyond the 
scope of the present paper. 



VI. CONCLUSIONS 

In conclusion, we analyzed the elementary excitations 
arising from the strong coupling of a photonic cavity 
mode with an intersubband transition of a single QW. 
We have shown that contrary to the current opinion 
Coulomb interactions can play crucial role in the sys- 
tem and lead to the qualitative changes of the nature of 
intersubband polaritons. We predict theoretically that 
strong coupling of the cavity mode occurs with collec- 
tive excitations, while single-particle excitations remain 
in the weak coupling regime. 

This work was supported by Rannis "Center of Ex- 
cellence in Polaritonics" and FP7 IRSES project "PO- 
LAPHEN". I.A.S. acknowledges the support from COST 
POLATOM program. O.K. acknowledges the help of 
Eimskip Foundation. 



Appendix A: Calculation of a polarization operator 
for non-interacting particles 

The generic form of the electron-hole polarization op- 
erator for non- interacting particles can be written as 



ino(a;,q) = 2 



dkdiy 



G(k + q,i^ + w)G(k,z/), (Al) 



7 



where G(k, v) denotes the Green function of particle with 
momentum k and energy v. q and uj correspond to the 
transferred momentum and energy, respectively. For the 
intersubband transition case the polarization bubble de- 
scribes the excitation process where electron is trans- 
ferred to the upper subband while the hole is created 
in the lower subband. Thus, the Eq. (jAl[) can be repre- 
sented as 



The integral (jA3p can be rewritten in the similar form 



no(cj,q) = -2 



2ir 
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2tt 



p(0) p(l) 
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2^ 



Bk + C + ij 
k 



k - ko + 



dk = (A8) 
dk = 
dk, 



Z7 



where E, 



(1) 



= A + te 



-k+q "'(k + q)V2m and E^^ = fi^k''/2m 

are dispersions of electrons in the excited and fundamen- 
tal subbands, respectively, A is an energy distance be- 
tween subbands and riok is the Fermi distribution in the 
fundamental subband which can be replaced by a step- 
like function at T = 0. Parameter 7 describes the lifetime 
of the excitation. 

First, let us rewrite Eq. (jA2p in the following form 

kdk 



2tt 



Ilo{uj,q) ^^^^2 hu- A- h^qy 2m- h-^kq COS 

(A3) 

where </> is an angle between two vectors k and q. The 
assumption that transferred momentum of photon q is 
small yields simple integration on (f> and fc, and the real 
part of polarization operator is 

5Rno(c.),_o = (A4) 

where n = ^ is 2D density of electron gas, A is QW 
area and we assumed 7—^0. 

Now we return to the case of non-negligible photon mo- 
mentum and calculate the imaginary part of polarization 
operator (jA3[) . The analysis of the denominator shows 
that integral on k in Eq. (|A3|) has poles only in certain 
angular range 



where B = lT?qcos(l)/m, C = hoj — A — lT?q^/2m, 
ko = C/B and the 0-dependence of B and kg should 
be taken into consideration. In the latest expression the 
new integration variable k = k — kp was used. Thus, 
comparing the Eqs. (jASp and (|A7[) . one can calculate 
numerically the imaginary part of polarization operator 
(|A3p with finite lifetime of excitations. One sees that 
integral (jASp can be separated into two parts. Detailed 
analysis shows that the first integral gives obtained pre- 
viously delta-function like absorption which docs not dc- 



7, while the second integral I Lor is proportional 
to 7 and has Lorentzian form. Therefore, for the sake of 
simplicity, in the purely analytical calculations it is pos- 
sible to use absorption spectrum as the sum of equation 
(jA6|) and phcnomenological Lorentzian broadening due 
to finite lifetime of excitation (r w /is) written in the 
formal 



nA 



hw- A 
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-jnA 



(fej- A)2+72- 



(A9) 



Finally, the imaginary part of polarization operator 
which describes optical absorption by intersubband tran- 
sition yields 



h?kpq/r 



< (h < arccos 



-{huj - A- Eq 



h?kpq/m 

(A5) 

Thereby, in this region one can find the imaginary part of 
polarization operator using residue theorem in the limit 
7^0: 



S>no(a;,g) = h{,oj,q) + lLor{(^) 



(AlO) 



The real part of polarization operator can be found by 
direct integration in Eq. 



sno(w,(7)~ 



Am 1 



{hw~A~Eq) tan( 



(A6) 

which describes a peak in the region where single-particle 
excitations exist, bounded by parabolas A -I- h^q^ /2m — 
h^kpq/m and A -I- h^q^ /2m + h^kpq/m. 

However, to describe the realistic absorption spectral 
function one needs to account for the finite non-radiative 
lifetime of intersubband excitations. It can be done using 
the Sokhatsky-Weierstrass theorem with non-zero 7 



mi,{u:,q) = ^^{{f^ - A - Eq)T 



T\ nnio-A-Eg)^ -AEpEq), 



(All) 



X + 



■dx 



Tr{x^ + 7^) 



dx 



(x2 4-72) 



dx. 



(A7) 



where " — " sign corresponds to the case hcu > A and 
" + " for hw <A. 

Consequently, the polarization operator Ho which 
stands for the simple intersubband bubble is described by 
the sum of real (Eq. (jAlip ) and imaginary (Eq. (|A10p ) 
parts. Finally, using it in the Dyson equations, the opti- 
cal response and elementary excitation spectrum can be 
found. 
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